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ABSTRACT 

Following its 2002 February eruption, V838 Mon developed a light echo that con- 
tinues to expand and evolve as light from the outburst scatters off progressively more 
distant circumstellar and/or interstellar material. Multi-filter images of the light echo, 
obtained with the South African Astronomical Observatory (SAAO) 1.0-m telescope 
between 2002 May and 2004 December, are analysed and made available electron- 
ically. The expansion of the light echo is measured from the images and the data 
compared with models for scattering by a thin sheet and a thin shell of dust. From 
these model results we infer that the dust is likely in the form of a thin sheet distant 
from the star, suggesting that the material is of interstellar origin, rather than being 
from earlier stages in the star's evolution. Although the fit is uncertain, we derive a 
stellar distance of ~ 9 kpc and a star-dust distance of ~ 5 pc, in good agreement 
with recent results reported from other methods. We also present JHKL and Cousins 
UBVRI photometry obtained at the SAAO during the star's second, third and fourth 
observing seasons post-outburst. These data show complex infrared colour behaviour 
while V838 Mon is slowly brightening in the optical. 

Key words: stars: individual: V838 Mon — stars: photometry — stars: circumstellar 
matter — stars: interstellar matter — stars: distances — stars: imaging 



1 INTRODUCTION 

V838 Mon appeared on the scene with a nova-like outburst 
of amplitude AV « 6 magnitudes in 2002 January (Brown 
2002) and its extraordinary subsequent evolution indicates 
that this object represents an entirely new type of variable 
(Munari et al. 2002, Evans et al. 2003). Besides its unique 
light curve, V838 Mon also underwent dramatic spectral 
evolution as the star rapidly expanded and cooled, trans- 
forming itself from an object vaguely resembling a mid-K 
giant into what Evans et al. (2003) propose to be the first 
known L supergiant. The emergence of a blue continuum in 
spectra obtained in 2002 October (Desidera & Munari 2002, 
Wagner & Starrfield 2002) suggests that V838 Mon may 
be a binary, although this early-type star could possibly be 
a line-of-sight companion. Various unconventional models, 
including merging main sequence stars (Soker & Tylenda 
2004) and an expanding red giant engulfing planets (Retter 
& Marom 2003), have been proposed to explain the star's 



behaviour. More recently, van Loon et al. (2004) argued that 
V838 Mon may be a single, low-mass AGB star experiencing 
thermal pulses while Munari et al. (private communication) 
and Tylenda, Soker & Szczerba (2004) both favour a young 
binary system with a massive progenitor for the outbursting 
component. These diverse interpretations of the 2002 phe- 
nomena indicate the depth of the mystery surrounding this 
unusual variable. 

Optical imaging 2 weeks after V838 Mon's 2002 Febru- 
ary eruption revealed the development of a light echo (Hen- 
den, Munari & Schwartz 2002); the result of light from the 
outburst being scattered into our line of sight by previously 
unseen circumstellar and/or interstellar material. We see the 
echo expanding as the light reaches progressively more dis- 
tant dust. The only other known Galactic light echo was 
associated with Nova Persei in 1901 (Couderc 1939) and 
hence the V838 Mon echo has received widespread atten- 
tion, particularly following the release of spectacular Hubble 
Space Telescope (HST) images^. 
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Aesthetics aside, the hght echo also provides impor- 
tant information about this enigmatic star. The observed 
expansion rate of the echo will eventually yield the nature 
of the scattering material, as well as the distance to the 
star and the distance between the star and the dust. Knowl- 
edge of V838 Mon's luminosity will allow stricter evaluation 
of the various models that attempt to account for the star's 
2002 eruptions. Furthermore, since the light echo illuminates 
progressively more distant material, each image provides a 
once-off map of the dust distribution on a given light echo 
paraboloid. Assembling these maps, taking into account the 
transformations required to correct for projection effects, 
would ultimately provide a 3-dimensional reconstruction of 
the scattering dust. 



2 OBSERVATIONS 

V838 Mon was included in the SAAO Infrared Service Pro- 
gramme following the initial outburst in 2002 January. Data 
from the first observing season were presented by Munari 
et al. (2002) and Crause et al. (2003). Subsequent JHKL 
measurements, also made with the IRP Mark II photometer 
on the 0.75-m telescope, are listed in Table 1. As for the 
first season data, a 36 arcsec aperture was used. Target ob- 
servations were interspersed with Carter infrared standards 
(Carter 1995) situated at similar airmasses and the data re- 
duced with SAAO software that interpolates between stan- 
dard stars to calibrate the sky-subtracted measurements. 
Given the infrared brightness of V838 Mon, contributions 
from faint field stars within the aperture are negligible. The 
typical uncertainty in the corrected magnitudes is 0.03 mag 
for the JHK data and 0.05 mag for the L band observations. 

All of the optical images and photometry presented here 
were obtained with the SAAO 1.0-m telescope and Ik x Ik 
STE4 CCD camera. The CCD was used without prebinning, 
resulting in a plate scale of 0.31 arcsec pixel"^. During 12 
observing runs between 2002 May and 2004 December, imag- 
ing was limited to moonless conditions with seeing typically 
better than 2 arcsec. The CCD frames were bias-subtracted, 
flat-fielded and trimmed using the ccdproc package within 
IRAF. 

Whenever possible, V838 Mon observations were brack- 
eted with red [{V — /) ~ 2] Landolt standard stars, but the 
extreme post-maximum colours [{V — i?) ~ 3 and {V — I) 
~ 6] required large extrapolations for the standard transfor- 
mations. PSF-modeling was performed with SAAO software 
and the resulting Cousins UBVRI photometry is listed in 
Table 2. The faintness of the star and the low sensitivity of 
the CCD in the U band made it impractical to obtain high 
signal-to-noise data and so U magnitudes are only quoted 
to one decimal place. Having obtained calibrated photom- 
etry at least once during most runs, it was possible to dif- 
ferentially correct additional observations made under non- 
photometric conditions using field stars located near V838 
Mon. These differential measurements are marked with stars 
in Table 2. Due to V838 Mon's extreme redness, the short 
exposures required to prevent saturation limited the accu- 
racy of the I band photometry as all comparison stars were 
severely underexposed. 

Throughout this paper we use Julian Dates in abbre- 
viated form, JD-2452000 (AJD). In Section 6 we introduce 



Table 1. SAAO JHKL photometry of V838 Mon obtained with 
the 0.75-m telescope between 2002 September 8 and 2004 De- 
cember 10. Dates are given as JD-2452000 (AJD) and horizontal 
bars separate the observations made during the second, third and 
early-fourth observing seasons. 
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529.666 
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6.95 


5.78 


5.09 


3.62 


933.627 


7.00 


5.82 


5.13 


3.68 


955.596 


7.03 


5.86 


5.16 


3.68 


969.598 


7.05 


5.86 


5.17 


3.73 


974.588 


7.08 


5.87 


5.17 


3.73 


1031.443 


7.10 


5.88 


5.22 


3.80 


1063.323 


7.08 


5.89 


5.25 


3.83 


1070.316 


7.13 


5.89 


5.24 


3.80 


1139.207 


7.20 


5.92 


5.27 


3.89 


1142.198 


7.20 


5.93 


5.27 


3.86 


1277.641 


7.47 


6.15 


5.45 


4.13 


1290.635 


7.43 


6.12 


5.43 


4.10 


1327.591 


7.41 


6.12 


5.44 


4.10 


1350.534 


7.49 


6.18 


5.47 


4.03 



day numbers that indicate the interval between the 2002 
February 5 (AJD 311) outburst peak and the observation in 
question. 



3 PHOTOMETRY 

V838 Mon's complete VRI and JHKL light curves are 
shown in Fig. 1. The U and B data are relatively sparse and 
hence were not included in this figure. Note that, to preserve 
detail in the JHKL curves, the vertical scale of the infrared 
(top) panel is twice that of the optical (bottom) panel. Arbi- 
trary alternating symbols were used to help distinguish the 
various bands. 

During the 4 month interval between the peak of the 
2002 outbursts and the star entering conjunction, V838 Mon 
faded by ~ 9, 8 and 4 magnitudes in V, R and / respectively. 
The infrared brightness peaked ~ 40 days after the opti- 
cal maximum and the trend of the fade amplitude decreas- 
ing with wavelength continued through the JHKL bands. 
Crause et al. (2003) interpreted this first season photomet- 
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Table 2. SAAO UBVRI photometry of V838 Mon obtained be- 
tween 2002 September 12 and 2004 December 12. Low signal-to- 
noise data are only quoted to 1 decimal place. Dates marked with 
a star indicate differentially corrected data from non-photometric 
nights and horizontal bars separate the observations made during 
the second, third and early-fourth observing seasons. 
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30 
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2.03 


5.72 


533.63* 
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-0.1 


2.74 
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ric behaviour as being indicative of a brief dust formation 
episode following the third optical outburst, however, this 
has not been confirmed spectroscopically. 

Although the star's brightness was decreasing in all 
bands in 2002 June (AJD ~ 430), V838 Mon reappeared 
in 2002 September (AJD ~ 530) with all but the V band 
showing a renewed increase. The R and I bands both bright- 
ened by ~ 0.8 magnitudes before approximately levelling off 
at 13.5 and 9.8 respectively while the JHKL Ught curves 
developed a distinct broad peak roughly centred on AJD 
600. The onset of this peak corresponds to the time when 
Evans et al. (2003) classified V838 Mon's spectrum as that 
of an L supergiant and the variation in the shape of this 
bump across the 4 infrared bands, becomming broader and 
less peaked from J through L, reflects the dominance of the 
cool component. 

Since the star remained near its minimum optical 
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Figure 1. Infrared (JHKL) and optical (VRI) light curves for 
V838 Mon, including first season data from Crause et al. (2003) 
and second, third and early-fourth season SAAO data reported 
in Tables 1 and 2. The magnitude scale for the infrared (upper) 
panel is twice that for the optical (lower) data. Note the bump 
at ~ 600 d visible in all bands except V , as well as the gradual 
brightening in the V band since AJD ~ 800. 



brightness throughout the second season, the constant pho- 
tometric contribution from the blue companion can be seen 
in the (U-B) column of Table 2. The steady brightening in 
the V band during the third season is a continuation of the 
behaviour reported by Crause (2003) when the star emerged 
from conjunction in 2003 September. A consequence of this 
brightening is that the {V — R) and {V — I) colours became 
bluer while the {B — V) colour continued to redden. 

The infrared (J—H), (H — K) and (K — L) colour curves 
for V838 Mon are shown in Fig. 2 and the various infrared 
colour-colour diagrams in Fig. 3. Asterisks mark data from 
the first season (reported by Munari et al. 2002 and Crause 
et al. 2003) while diamond symbols indicate those presented 
in this paper. Data from the second, third and early-fourth 
seasons are shown as black, grey and open symbols respec- 
tively. The gaps roughly centred on AJD 460, 830 and 1200 
(2002, 2003 and 2004 July) are due to the star being in con- 
junction with the Sun. 

Figs 2a and 2b show the star recovering from its reddest 
(J — H) and (H — K) colours during the second season (2002 
September to 2003 May; AJD 526-776) while Fig. 2c reveals 
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Figure 2. (J — H), {H — K) and {K — L) colour curves for V838 Mon. Asterisks indicate published data from the first observing season 
while black, grey and open diamonds mark new data from the second, third and early-fourth seasons respectively. 
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Figure 3. (J — H) vs {H — K), (J — H) vs (K — L) and {H — K) vs {K — L) colour-colour diagrams showing the complex infrared 
evolution of V838 Mon. Symbols as in Fig. 2. The locus of blackbody colours reddened by Ay = 2.8 mag, as well as curves illustrating 
the temporal evolution of the data, are included in each panel. See Crause et al. (2003) for details of the 2002 infrared behaviour. 



that the {K — L) colour only reached its reddest value about 
170 d after the (J - H) and {H - K) colours did. Since the 
third season (2003 September; AJD 895), the {J - H) curve 
showed a redward trend while the {K — L) colour became 
slightly bluer. 

V838 Mon's complex colour behaviour is further illus- 
trated in Fig. 3. Schematic curves in the lower-right corners 
of the 3 panels indicate the temporal evolution of the colours 
in each plane. The first and early-second season photometry 
shows trends in all three colour-colour planes that are gen- 
erally consistent with the star's post-outburst evolution to- 
wards later spectral types and therefore lower temperatures. 
The decrease in temperature is illustrated here by compari- 
son with the locus of blackbody colours, after adjustment for 
reddening at the level of Ay = 2.8 magnitudes found to be 
appropriate for V838 Mon by Munari, Desidera & Henden 
(2002). Although blackbody colours are only approximate 
to the colours of stellar atmospheres, the early evolution 
of V838 Mon is seen to run roughly parallel to the black- 



body lines. The comparison further highlights the shift to- 
wards higher {K — L) colours in the middle first season data 
noted by Crause et al. (2003), and interpreted by them as 
evidence for a short-lived dust formation event. Thereafter, 
throughout the remainder of the second and third seasons, 
and now with the availability of early-fourth season pho- 
tometry, we see the star following complex paths and possi- 
bly the development of loops in the colour-colour diagrams. 
We would hope to elucidate these variations once late-event 
spectroscopy of the object becomes more readily available 
in the literature, e.g. the spectrum obtained by Evans et al. 
(2003), announcing the star as an L-type supergiant, was 
obtained early in the second season of observations (2002 
October 29; AJD 577). Our nearest photometric observation 
(at AJD 574.6; see Table 1) records the colours of the star 
to be [J -H) = 1.24, {H - K) = 0.99 and [K - L) = 1.38. 
In Figure 3, this places the star beyond the position of red- 
dest (J — H) colour which likely occurred when V838 was 
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Figure 4. The spectral energy distribution of V838 Mon was 
calculated at various epochs. Since the second observing season, 
the star has faded in the near infrared while brightening in the 
optical. The AJD numbers shown are average values for the com- 
bination of optical and infrared measurements used to produce 
each curve. 



in conjunction with the Sun in mid-2002, and during a time 
of decreasing (J — H) and {H — K) colour. 

The brightening seen in the V band during the third and 
fourth seasons is in contrast to the near infrared magnitudes 
that were fading during the same time period; see Figure 1. 
The spectral energy distribution of the star was calculated at 
various epochs, corresponding to 2002 September 12 (AJD 
530), 2003 April 20 (AJD 750), 2004 May 12 (AJD 1138) 
and 2004 December 10 (AJD 1350). These curves are shown 
in Fig. 4 and illustrate changes in the flux distribution be- 
tween the time of the faintest V^-band measurement (AJD 
530) and our most-recent observations centred at AJD 1350. 
In particular note the « 0.4 dex (factor of 2.5 in flux, or ap- 
proximately 1 magnitude) decrease in flux across the near 
infrared bands, and the nearly similar increase in flux seen 
in the B and V bands. 



4 LIGHT ECHO IMAGES 

The SAAO 1.0-m telescope was used to collect BVRI im- 
ages of the light echo between 2002 May and 2004 Decem- 
ber. Median averages of 3 or more frames were produced 
whenever possible to eliminate cosmic ray strikes and to 
increase the signal-to-noise ratio. When only 1 image was 
available, the IRAF task imedit was used to manually re- 
move cosmic rays and, when necessary, images were ad- 
justed with rotate to compensate for any differences in 
the camera orientation from run to run. Individual images 
were finally cropped from 5.3 x 5.3 arcmin to 2.5 x 2.5 
arcmin. Fig. 5 shows representative BVRI images chosen 
to illustrate the yearly expansion of the light echo. The 
2002-2004 sequences for each filter were scaled for expo- 
sure times and then displayed on the same intensity scale 
to facilitate direct comparison. The BVRI frames obtained 
under the best seeing and transparency conditions for each 
epoch are available in FITS format and may be downloaded 
from ftp:/ /ftp.saao.ac.za/pub/lisa/V838images/ Dates, ex- 



posure times and seeing details are given in Table 3 and in 
a 'readme' file on the ftp site. 

The standard model for the light echo phenomenon in- 
volves a short duration light pulse that propagates outwards 
and scatters off surrounding dust. In the case of V838 Mon 
however, the star experienced an optical maximum that 
lasted about 70 d (see Fig. 1) and during this time the 
light source also changed significantly as the star rapidly ex- 
panded and cooled. Since V838 Mon was at its hottest dur- 
ing the 2002 February outburst, the light echo was brightest 
in the U and B bands when it was discovered (Henden et al. 
2002). The brightness of the echo decreases due to geomet- 
ric dilution and the colour evolution follows that of the star, 
while also being influenced by the scattering process that is 
dependent on the dust properties. 

Fig. 5 shows that the echo faded more rapidly at shorter 
wavelengths than in the R and / bands and so the dominant 
factor in the echo colour evolution appears to be the be- 
haviour of the light source. For particles much smaller than 
the wavelength of light (the Rayleigh regime), the scattering 
efficiency is highly wavelength dependent - proportional to 
A^**. If this were the case, we might expect the light echo to 
have been much fainter and to have faded faster in the red. 
This argument suggests that the dust particles producing 
the echo are within the Mie scattering regime (i.e. the par- 
ticle sizes are comparable to the wavelength of light), which 
is consistent with the properties of interstellar dust grains 
(Glass 1999). 



5 ESTIMATING THE DISTANCE FROM THE 
LIGHT ECHO 

The geometry of a light echo dictates that, provided mate- 
rial remains present to scatter the light, the echo will appear 
circular regardless of the physical distribution of the scatter- 
ing material; see Tylenda (2004) for further discussion. Thus 
to measure the size of the light echo we fitted a circle to the 
outermost edge of the refiection nebula at each epoch and 
used the radius as a measure of the echo's angular size. We 
measured our V band images as they usually had the high- 
est signal-to-noise ratio in the nebula. The data presented in 
Table 3 show that the V band images had a mean seeing of 
~ 1.6 arcsec, with a maximum difference of ~ 1 arcsec across 
the various epochs. These variations in seeing contribute to 
the overall uncertainty of the circular fit, which is domi- 
nated by the inherent difficulty in establishing the 'edge' of 
the light echo. While the nebula is sharply defined in the 
early images, it becomes more diffuse at later epochs and 
hence more difficult to measure. Accordingly, we estimate 
the uncertainty in the measured radii to be ~ 3 percent, 
which gives an absolute error in the radii that scales pro- 
portionally to the light echo radius as it expands with time; 
see the error bars in Fig. 6. Table 4 gives the measured light 
echo radius at each epoch, along with the day numbers that 
denote the interval between each observation and the time 
of the star's maximum brightness (2002 February 5; AJD 
311), when the light echo is believed to have been triggered. 

In measuring the radii, we found that the centre of the 
light echo systematically shifted away from V838 Mon, to- 
wards the NE (in the upper-left direction in Fig. 5). This 
shift tends to make the echo seem asymmetric, but merely 
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Figure 5. B, V, i? and / images of tlie expanding light eclio obtained with the SAAO CCD on 2002 May 1 (AJD 396), 2003 May 21 
(AJD 781) and 2004 May 7 (AJD 1133). Differences in the appearance of stars in the 3 images for each band are due to seeing variations 
while the colour evolution of the light echo clearly follows that of V838 Mon. North is up and east to the left in each of the 2.5 X 2.5 
arcmin frames. The measured light echo diameter is indicated by a horizontal bar in each of the V images. See Table 3 for further details. 



ofT-setting the circle being fitted compensates for this effect. 
Such a shift suggests that the scattering dust is in the form 
of a sheet that is inclined relative to the line-of-sight; the 
measured Right Ascension and Declination offsets in units 
of arcsec are listed in Table 4. These offset values have un- 
certainties of a few arcsec, but clearly follow a linear trend 
that allows us to derive the inclination angle of the implied 



sheet once an estimate of the distance to the star has been 
obtained (for details see Tylenda 2004). 

Tylenda (2004) modeled the light echo evolution for 
dust distributed in thin interstellar sheets and in thin cir- 
cumstellar shells. Both cases initially produce an expanding 
echo that cannot be distinguished by model fits to the ob- 
servations. However, the models predict a divergence with 
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Table 3. Dates, exposure times (seconds) and seeing measurements (arcsec) for the best BVRI images for each epoch. These were 
obtained with the SAAO 1.0-m telescope between 2002 May 1 and 2004 December 11. During an observing run, images for the various 
filters were often obtained on different nights; the UT dates given below are for the V images. The exact dates for all the frames in a 
given set differ by up to 3 days, details may be found in the headers of the files available from ftp: / /ftp.saao.ac.za/pub/lisa/V838imagcs/_ 



UT Date 


AJD 


B band 


B band 


V band 


V band 


R band 


R band 


I band 


/ band 






exposure 


seeing 


exposure 


seeing 


exposure 


seeing 


exposure 


seeing 


2002 May 01 


396 


1 X 600 


1.4 


1 X 300 


1.2 


1 X 20 


1.3 


1x2 


0.8 


2002 Jun 02 


428 






1 X 300 


2.1 


1 X 100 


2.2 


1 X 30 


2.0 


2002 Sep 16 


534 


3 X 300 


2.3 


3 X 600 


2.0 


3 X 300 


1.8 


3 X 60 


1.7 


2002 Oct 31 


579 


1 X 600 


2.3 


1 X 600 


1.9 


3 X 300 


1.8 


5 X 20 


1.4 


2003 Jan 26 


666 


3 X 1200 


1.5 


3 X 900 


1.0 


3 X 300 


1.5 


5 X 10 


1.1 


2003 Apr 19 


749 


3 X 900 


1.9 


3 X 900 


1.8 


3 X 300 


1.4 


9 X 15 


1.5 


2003 May 21 


781 


3 X 1200 


1.5 


2 X 600 


1.3 


3 X 300 


1.1 


9 X 10 


0.9 


2003 Sep 12 


895 


3 X 900 


2.3 


3 X 600 


1.6 


3 X 250 


1.4 


5 X 10 


1.0 


2003 Nov 03 


947 






3 X 900 


1.5 


3 X 300 


1.4 


7 X 20 


1.2 


2004 Mar 24 


1089 


1 X 600 


1.7 


1 X 200 


1.5 


1 X 200 


1.4 


1 X 15 


1.3 


2004 May 07 


1133 


1 X 1800 


2.0 


3 X 900 


1.4 


3 X 300 


1.3 


5 X 10 


0.9 


2004 Dec 11 


1351 


3 X 1800 


2.0 


5 X 1200 


1.3 


3 X 600 


1.6 


11 X 25 


1.6 



time - shells of circumstellar dust will eventually result in 
the light echo contracting and disappearing, while expansion 
due to interstellar dust will decelerate but still continue as 
the echo gradually fades. Resolving the nature of the scat- 
tering material justifies the continued interest in imaging of 
the nebula. For both the dust distribution models consid- 
ered here, the angular expansion of the light echo can, in 
principle, be used to determine the stellar distance and the 
distance of the scattering material from the star. 

Fig. 6 shows our measured radii versus the time since 
the peak of the 2002 February outburst (AJD-311). To these 
data we fitted models for scattering by a single thin sheet 
(+ sign) and a single thin shell (— sign) of material, making 
use of the equation 9 = 1/d \/2rct ± c^f^, where 9 is the 
angular radius of the echo, d is the distance to the star, r is 
the distance between the star and the dust, c is the speed 
of light and t is the time since the peak of the outburst. 
This form is equivalent to that given by van Loon et al. 
(2004), other than invoking the small angle approximation 
as the angles involved will always be sufficiently small that 
tan 9x9. 

minimisation yielded r and d values of 4.8 ± 1.8 pc 
and 8.9 ± 1.6 kpc for the sheet model and 9.4 ± 9.8 pc and 
12.0 ± 6.5 kpc for the shell model, respectively. The model- 
derived results are sensitive to uncertainties, or errors, in 
the measured radii^, with the shell case being additionally 



complicated by the need to provide a reasonable initial value 
for the radius r so as to avoid obtaining a negative term 
within the square root term of the expansion equation. 

Tylenda (2004) discusses the intrinsic uncertainty of the 
fits due to the shallow, extended nature of the minimum, 
but it is encouraging to note that both models produced val- 
ues for the distance to V838 Mon that are consistent with 
lower limits proposed by various authors (e.g. Kipper et al. 
2004, Wisnewski et al. 2003, Bond et al. 2003) and in agree- 
ment with values for the distance of 9.2 kpc and 8 ± 2 kpc 
derived recently by Lynch et al. (2004) and Tylenda (2004), 
respectively. Fig. 6 suggests that the expansion of the echo 
may be entering a regime where it has become possible to 
distinguish between the two dust distribution models owing 
to the significant divergence in the model-predicted light 
echo radii. 

Having obtained a distance d — 8.9 ± 1.6 kpc from 
the sheet model, we estimated the inclination angle (a) 
of the dust sheet's normal to the line-of-sight. Following 
Tylenda (2004), we applied a linear fit to the measured 
RA-offset of the echo-centre versus time (since outburst) 
and derived a slope of = LOl x 10"^ arcsec day~^. Since 
tan a = 9cd/ct, appropriate unit conversions (arcsec day~^ 
into radians day ~^ and then radians into degrees) yielded 
a 27°. 



^ van Loon et al. (2004) measured the diameter of the light echo 
in images obtained with various telescopes, including the HST 
and the SAAO 1.0-m, and attempted to fit thin sheet and thin 
shell models to these data. As they were unable to find solutions 
for either model that fitted both the early and later epochs, they 
relaxed the assumption of a thin scattering medium and derived 
a lower limit of > 5.5 kpc for the distance to V838 Mon. However, 
comparison between the diameters given by van Loon et al. (see 
their table 1) and the radii listed by Tylenda (2004; see table 1 
and 2) and in Table 4 of this paper reveals that the van Loon et al. 
diameters are in error, underestimating the actual echo diameters 
by a factor of ~ 2.5. We also note that the UT dates estimated 
for their third and sixth epochs (SAAO images) are in error by 
16 and 1 d respectively. 



6 DISCUSSION 

Lynch et al. (2004) model 0.8 - 2.5, 2.1 - 4.6 and 3 - 14/im 
infrared spectra with a central star and a 2-component ex- 
panding circumstellar shell. They treat H2O as the dominant 
molecular absorber and emitter in the envelope which atten- 
uates stellar emission in the near infrared and emits ther- 
mally in the infrared. H2O and CO transmission through 
the shell dominates in the 2.5 — 4.5/im region and the pho- 
tosphere makes a minimal contribution beyond 4.5/im. 

Their model does not include silicate dust and they 
point out that the 2002 February/March extinction can be 
explained with a combination of molecular Rayleigh scat- 
tering and molecule formation. However, as early as 2002 
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Table 4. Measured light echo radii for each of the V band im- 
ages and their respective dates, as well as a time representing 
the number of days since the peak of the 2002 February outburst 
(AJD-311). Right Ascension («) and Declination (5) offsets be- 
tween V838 Mon and the centre of the echo arc also given. 



UT Date 


AJD 


Time 
(d) 


Radius 
(arcsec) 


a-offset 
(arcsec) 


(5-offset 
(arcsec) 


2002 May 01 


396 


85 


19.5 


-1 


1 


2002 Jun 02 


428 


117 


21.4 


-0 





2002 Sep 16 


534 


223 


31.9 


-2 





2002 Oct 31 


579 


268 


34.1 


-2 


2 


2003 Jan 26 


666 


355 


39.4 


-3 


2 


2003 Apr 19 


749 


438 


44.3 


-4 


3 


2003 May 21 


781 


470 


46.5 


-5 


3 


2003 Sep 12 


895 


584 


51.8 


-6 


3 


2003 Nov 03 


947 


636 


53.9 


-7 


3 


2004 Mar 24 


1089 


778 


58.6 


-7 


5 


2004 May 07 


1133 


822 


61.4 


-8 


7 


2004 Dec 11 


1351 


1040 


70.1 
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Figure 6. Measured light echo radius (arcsec) versus the time 
(d) since the peak outburst in 2002 February. Best-fit models for 
scattering material distributed in a thin sheet and a thin shell 
are shown and the values for the stellar distance, d, and star-dust 
separation, r, derived from each fit are listed. Error bars indicate 
the estimated 3 percent uncertainty in the measured echo radii. 



January, their spectra showed a slightly elevated 10/xm re- 
gion and by 2003 January they report that the 8 — 13/im 
flux had increased 10 times more than that in the near in- 
frared. While the development of this infrared excess could 
be due to gas-phase molecules in the circumstellar envelope, 
a large amount of dust-precursor material was present in 
the form of molecular SiO and in 2003 February a fairly 
strong feature near lO/im was tentatively identified as sil- 
icate emission with a central absorption feature (Sitko et 
al. 2003). The timing of this observation roughly coincides 
with the end of the broad peak noted in the infrared light 
curves (Fig. 1) and the time of the maximum {K — L) colour 
of V838 Mon (Fig. 2c). This confirms that at least some of 
the infrared photometric activity was caused by the produc- 



tion of dust, although the extent to which dust formation 
occurred throughout the outburst remains in question. 

Light echoes have been detected in several extra- 
galactic supernovae (Chevalier 1986) , the best studied being 
that of SN 1987A (Xu, Crotts & Kunkel 1995). However, 
prior to V838 Mon, Nova Persei 1901 (Couderc 1939) was 
the only Galactic variable known to have developed a light 
echo. 

While most observed light echoes have been due to in- 
terstellar dust, it is conceivable that the scattering could 
take place within a circumstellar envelope. The evolution- 
ary status of V838 Mon remains uncertain, but circumstellar 
dust would be expected if the star has undergone multiple 
outbursts in the past. The discovery of extended cool shells 
surrounding V838 Mon reported by van Loon et al. (2004) 
and detections by both 2MASS and IRAS (Kato 2002) may 
imply that the progenitor was an evolved object. 

However, if V838 Mon is indeed kinematically associ- 
ated with a B3V companion (Desidera & Munari 2002, Wag- 
ner & Starrfield 2002), the system has to be young, indi- 
cating a massive progenitor for the outbursting component. 
This outcome is consistent with Tylenda et al. (2004) who 
claim the progenitor was most likely a 5 — 10 M© dwarf, and 
who also reject the link between V838 Mon and the extended 
cool emission found by van Loon et al. (2004). Reddening 
of the B3V star indicates that the system is ~ 10.5 kpc 
distant (Munari et al. 2002) and although the star is situ- 
ated towards the Galactic anti-centre, its Galactic latitude 
of 1.05° places it within 200 pc of the plane, consistent with 
the observed extinction being due to interstellar material. 

Also, Wisnewski et al. (2003) and Desidera et al. (2004) 
conclude that the polarization observed in V838 Mon is 
mostly of interstellar origin. A small intrinsic component was 
present during the 2002 February maximum, but decreased 
rapidly after that. Furthermore, deep polarimetric imaging 
of the light echo (Desidera et al. 2004) shows a dipole struc- 
ture that would not be consistent with scattering from a 
homogeneous circumstellar medium. 

While the real dust distribution that results in the light 
echo is likely to be complex. Fig. 6 shows that the observed 
light echo expansion can be reproduced by simple models 
that assume a thin sheet or a thin shell of scattering mate- 
rial. Tylenda (2004) analysed the HST images of the fight 
echo and found that a sheet of material matched the mea- 
sured radii better than a shell model could. Only a shell of 
implausibly large radius, hence approximating a flat sheet of 
material, could be made to fit the data. Our larger fit-errors 
for the shell (105 and 54 percent for r and d respectively, 
compared to 37 and 18 percent for the sheet) also indicate 
that the sheet model more successfully represents the data. 

We expect the nature of the scattering material will 
be established by the end of the fourth observing season as 
the sheet and shell models will have separated significantly 
by that time and hence allow discrimination between the 2 
possibilities. With the echo fading due to geometric dilution, 
high signal-to-noise images will be required to avoid under- 
estimating the radius as that would bias fits towards the 
shell model. 

Future measurements of the radius will also help to re- 
fine the distance as later points will improve the model fit 
and hence provide better estimates of d and r. At a much 
later stage, the observed shift between the echo centre and 



V838 Mon: light echo evolution and distance estimate 9 



the star could be used, in conjunction with the radius mea- 
surements, as an alternative method to determine the star's 
distance (Tylenda 2004). 



7 SUMMARY 

V838 Mon seemed to have settled down after the dramatic 
2002 events, but the photometry presented in Tables 1 and 2 
indicate that the properties of the object continued to evolve 
(see Fig. 1). During the second season, the star remained 
faint in the V band but brightened by ~ 0.8 mag in R and 

1 and developed a broad, flat bump in the JHKL bands 
that peaked around AJD 600. The JKL infrared magnitudes 
then faded slowly during the third season while the H and 
/ bands remained fairly constant and the V band began 
brightening. The complex infrared evolution shown in Figs 

2 and 3 and changes in the overall flux distribution shown 
in Fig. 4 appear consistent with a dust dispersal scenario. 

Fig. 5 consists of sample BVRI light echo images span- 
ning a 2 year interval; see Table 3 for all image details. Our 
current set of images may be downloaded from the SAAO 
ftp site; see Section 4 for details. 

We measured the radius of the expanding light echo 
in our V band images and fitted models for a thin sheet 
and a thin shell of scattering material. While the thin sheet 
model gave a somewhat better fit to the light echo radii, 
both models yield distances (8.9 ± 1.6 kpc and 12.0 ± 6.5 
kpc for the sheet and shell respectively) that are consistent 
with estimates from other methods - see Fig. 6 and Table 4. 
However, the colour evolution of the light echo, the higher 
quality of the model fit for a sheet of dust, the observed 
shift in the echo centre relative to the star, results from 
polarimetric studies and the position of the variable only 
slightly above the Galactic plane all suggest that the dust 
illuminated by the outburst of V838 Mon is more likely to 
be of interstellar origin. 

Establishing the nature of the dust may provide evo- 
lutionary information while the distance is essential for ac- 
curately determining the absolute luminosity of this myste- 
rious variable. Continued imaging thus remains a priority; 
with time the dust distribution will be revealed and the bet- 
ter constrained fit will improve the distance estimate. 
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